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Classification of complex wave functions of infinite variables is an important problem since it is 
related to the classification of possible quantum states of matter. In this paper, we propose a way 
to classify symmetric polynomials of infinite variables using the pattern of zeros of the polynomials. 
Such a classification leads to a construction of a class of simple non-Abelian quantum Hall states 
which are closely related to parafermion conformal field theories. 
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I. INTRODUCTION 
A. Functions of infinite variables 

One of the most important problems in condensed mat- 
ter physics is to understand how particles are organized 
in a ground state. Almost all the low energy properties 
of the system are determined by such an organization. 
Mathematically, the ground state of N particles is de- 
scribed by a wave function - a complex function of N 
variables <i>(ri,-- - ,rjv), where is the coordinate of 
the i**^ particle. Thus, the problem to understand the 
patterns of the many-particle organizations (or in phys- 
ical terms, to understand the quantum phases of many- 
particle systems) is to classify the complex wave functions 
<I>(ri, • • • ,rjv) in the N ^ oo limit. 

Such a classification problem is one of most fundamen- 
tal problems in physics, since it determines the possible 
quantum phases of many-particle systems. Due to the 
success of Landau symmetry breaking theory in describ- 
ing phases and phase transitions,^ for a long time, physi- 
cists believe that phases of matter are classified by their 
symmetry properties. Mathematically, this is equivalent 
to believing that the wave functions are classified by their 
symmetry properties, for example whether the wave func- 
tion is invariant under translation $(ri) — > '^{vi + a) or 
not. Under such a belief, the wave functions with the 
same symmetries are grouped into one class and such a 
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class represents a single phase of matter. This is why 
group theory becomes an important mathematical foun- 
dation in physics. 

However, after the discovery of fractional quantum 
Hall (FQH) states, ^'•^ it was realized that symmetry is not 
enough to classify all the possible organizations encoded 
in the wave functions $(ri, • • • , rjv)- This is because the 
wave functions that describe different FQH states have 
exactly the same symmetry. Thus the wave functions of 
FQH states contain new kinds of organizations of parti- 
cles that has nothing to do with symmetry.^'^ The new 
organizations of the particles are called topological or- 
ders. 

Intuitively, what is new in the FQH wave functions 

is that the wave functions contain a long-range quantum 
entanglement.^'^ This is why the FQH wave functions de- 
scribe new states of matter that cannot be described by 
symmetries. The wave functions with long-range entan- 
glements and the corresponding topological orders not 
only appear in FQH systems, they also appear in various 
quantum spin systems. Understanding this new class of 
wave functions and the resulting new states of matter is 
currently a very active research direction in condensed 
matter physics. 

To gain a deeper and more precise understanding of 
topological orders and the associated long-range entan- 
glements, we need to solve the related mathematical 
problem of classifying $(ri,--- ,rAr) in the N ^ oo 
limit. This is a difficult problem which is not well studied 
in mathematics. The problem is not even well defined. 
But this does not mean that the problem is not impor- 
tant. It is common not to have a well defined problem 
when we wander into an unknown territory. The first 
task of knowing the unknown is usually to come up with 
a proper definition of the problem. 

In this paper, we will not attempt to classify generic 
complex wave functions <I>(ri,-- - jV^). We limit our- 
selves to a simpler problem of trying to classify FQH 
states and their topological orders. (For a review on 
topological order in FQH states, see Ref. 29,30.) The 
corresponding mathematical problem is to classify sym- 
metric and anti-symmetry polynomials of N variables 
^{zi,--- ,Zn) in the N ^ oo limit. We will first try 
to come up with a physically meaningful and mathemat- 
ically rigorous definition of the problem. Then we will 
solve the problem in some simple cases. This leads to 
a class of "simple" (anti-)symmetric polynomials which 
corresponds to a class of "simple" FQH states. The 
constructed FQH states include both Abelian and non- 
Abelian FQH states. "-^^ 



II. FQH STATES AND POLYNOMIALS 

A. FQH wave functions 

First we would like to give a brief review on FQH the- 
ory. A FQH state is a quantum ground state of two di- 



mensional electrons in a magnetic field. Such a quantum 
state is described by a complex wave function 

^ixi,yi,X2,y2, ■ ■ ■ ,xn,Vn) 

where [xi, yi) are the coordinates of the i^^ electron and 
N the total number of electrons. In the strong magnetic 
field limit, if the filling fraction z/ of a FQH state is less 
than 1, then all the electrons are in the lowest Landau 
level. In this case, the ground state wave function has 
the following form 

* = ^{zi, ■ ■ , ZN)e~i ^•'^i l""*!^ 

where Zi = Xi + iyi and ^{zi, - ■ • , zn) is a holomorphic 
function of Zi [ie ^{zi, ■ ■ ■ ,Zm) does not depend on z*]. 
Since ^{zi, ■ • ■ , ^^v) has no poles, • • ■ , zn) is a poly- 

nomial of Zi's. 

Due to the Fermi statistics of the electrons, 
^{z\, • • • , Zn) must be an antisymmetric polynomial. If 
we assume the electrons to have a Bose statistics, then 
$(2;i,--- ,Zn) must be a symmetric polynomial. So to 
understand the phases of FQH systems is to classify an- 
tisymmetric or symmetric polynomials. 

It turns out that for every antisymmetric polynomials 
^aiiti-sys('2i, ' ' ' ,Zn), One Can uniquely construct a sym- 
metric polynomials #sys(^i,--- .zn) = ^'t^'^fc-'.;)"^ " 
Thus classifying anti-symmetric polynomials and classi- 
fying symmetric polynomials are almost identical prob- 
lems. So in this paper, we will assume electrons to have a 
Bose statistics and concentrate on classifying symmetric 
polynomials. 

For a system of N bosonic electrons, which symmetric 
polynomial will represent the ground state of the system? 
It will depend on the interaction between the electrons. 
If the interaction potential between two electrons has a 
(5-function form 

Fi(2i,Z2)='5(zi-Z2) (1) 

then the ground state is described by the symmetric poly- 
nomial 

^1/2 = J{{Zi - Zjf. 
i<j 

Such a state has a vanishing total potential energy Vtot = 
0, where 

Vtot = J l[d^Zi^*{{Zi})J2ViZi,Zj)^{{Zi})- 

i i<j 

The vanishing of the total potential energy Vtot requires 
that the wave function '^{{zi} to be zero as zi — > Z2- 
Since the average energy Vtot > for any wave func- 
tions, the vanishing Vtot for $1/2 indicates that $1/2 is 
the ground state. 

If the interaction potential between two electrons is 
given by^^ 

V2(^i, Z2) = vod{zi - Z2) + V2d^,.S{zi - Z2)dl (2) 
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with t;o > and V2 > 0, then the ground state will be 

^1/4 = n(^' - ^j)''- 

i<j 

For interaction (2), the vanishing of the total poten- 
tial energy Vtot not only requires that the wave function 
^'({Zj} to be zero as Zi ^ Z2, it also requires ^'({2;,} to 
vanish faster than (zi — 2:2)^ as zi — > Z2- This means 
that the symmetric polynomial must have a fourth or- 
der zero as Zi — » Zj. One such polynomial is given by 
^1/4 = ni<j(-2^i ^j)^ which has the lowest total power 
of Zi's. 

More complicated ground states can be obtained 
through more complicated interactions. For example, 
consider the following three-body interaction between 
electrons^'''^* 

Vpf{zi,Z2:Zs) = (3) 

S[vo6{zi - Z2)6{z2 - Z3) - viS{zi - Z2)d^*S{z2 - 23)5^3] 

where S is the total symmctrization operator between 
Zi, Z2, and Z3. Such an interaction selects the following 
symmetric polynomial^^ 

$P/ = a( — ) TT(zi - Zj) 

\Zi- Z2 Z3 — Zi Zn-1 - Zn^ ^J-. 

to describe the ground state (which has a vanishing to- 
tal potential energy Vtot)- Here A is the total antisym- 
metrization operator between zi, - ■ ■ ,zn. 

Three symmetric polynomials ^1/2, ^i/i, and i>p/ 
contain different topological orders and correspond to 
three different phases of an iV-electron system in the 
N 00 limit. They are the filling fraction v = 1/2 
Laughlin state, ^ the filling fraction = 1/4 Laughlin 
state, and the filling fraction v — 1 Pffafian state. We 
would like to find a classification of symmetric polyno- 
mials such that the above three symmetric polynomials 
belong to three different classes. 

B. Ideal Hamiltonian and zero-energy state 

The above three examples share some common prop- 
erties. The Hamiltonians described by the interaction 
potentials Vi, V2, and Vpf arc all positive definite and 
contain zero-energy eigcnstates. The zero-energy eigen- 
states of Vi, V2, and Vpf are known and are given by 
^1/2) ^1/4) and $p/. Thus $1/2, $1/4, and $p/ are 
exact ground states of the corresponding Hamiltonians. 
Since the interaction potentials are constructed from 5- 
functions and their derivatives, the exact ground states 
(the zero-energy states) for such type of potentials are 
characterized by the pattern of zeros, ie the orders of ze- 
ros of the ground state wave function as we bring two or 
more electrons together. 

In this paper, we will concentrate on such ideal Hamil- 
tonians and their exact zero-energy ground states. From 



this point of view, classifying FQH states corresponds to 
classifying patterns of zeros in symmetric polynomials. 
In other words, for each pattern of zeros, we can define 
an ideal Hamiltonian such that the symmetric polynomi- 
als with the given pattern of zeros will be the zero energy 
ground states of the Hamiltonian. Such symmetric poly- 
nomials will describe a phase of a FQH system, provided 
that the Hamiltonian has a finite energy gap. 

Clearly, for the ideal Hamiltonians, apart from the 
zero-energy ground states, other eigenstates of the Hamil- 
tonian all have non-zero and positive energies. But this 
does not imply the Hamiltonian to have a finite energy 
gap. Only when the minimal energy of the excitations has 
a finite non-zero limit when electron number approaches 
to infinity, does the Hamiltonian have a finite energy gap. 
So to classify the FQH states, we not only need to classify 
the patterns of zeros and the associated ideal Hamiltoni- 
ans, we also need to judge if the constructed ideal Hamil- 
tonian has a finite energy gap or not. At the moment, 
wc do not have a good way to make such a judgement. 
So here, we will concentrate on classifying the patterns 
of zeros and the associated zero energy states. 

III. PATTERN OF ZEROS 

A. Derived polynomials and their Dab 
chciracterization 

In order to classify translation invariant symmetric 
polynomials of N variables ^(zi, • • • , zat) in the N ^ 00 
limit, we need to define the polynomials for any N. The 
key in our definition is to introduce "local conditions". 
Those local conditions apply to polynomials of any num- 
bers of variables. 

Prom the discussion in the Section II A, we see that 
one way to implement the local condition is to let one 
variable to approach another and to specify the power of 
the zero 

$(21,- •• ,Zn)\zi^Z2 ~ (21-22)^". 

We would like to stress that the above local condition is 
consistent with the translation invariance of the polyno- 
mial. 

However, Dn does not contain all the information 

that is needed to specify various interesting polynomials. 
To implement more general local conditions, we need to 
bring three or more variables together and specify the 
patterns of zeros. ^^'^^ 

To describe the patterns of zeros in a systematic way, 
we obtain from a polynomial $(21, • • • , Zn) another poly- 
nomial P' by letting zi — > 22: 

^{Zl, ■■ ■ , Ziv)Li^Z2=z(2) 

~(21-Z2)^"P'(Z(2\Z3,-- - ,Zn) 

Here ^ means equal up to a non-zero complex constant. 
The value of Dn encodes a part of the local conditions. 
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Then we let Z3 



, Zn)- 



P'{z^^\zs, ■■■ ,ZN)r^{zs- z^'Y''P"{z^'\z,, zn), 

where z^^^ = z^'^\ In this way, we obtain a new polyno- 
mial P"(z("^\ 2:4, •• • , Zn)- In general wc obtain P({z|°''}) 
where z^"-^ is a type-a variable obtained by fusing a Zi- 
variables together. Note that z^^^ = Zi is the original 
variable. If we view Zi = z'l^^ as coordinates of electrons, 
then z|"'' arc coordinates of bound states of a electrons. 
We will call such a bound state a type-a particle. 

P{{zl°'^}) is a symmetric polynomial that is symmetric 
between variables of the same type. It satisfies certain 
local conditions and form a Hilbert space. The polyno- 
mial P({z-"^}) is also called a derived polynomial since 
it is obtained from <I>({zj}) by fusing variables together. 

The general local conditions on ^{{zi}) are specified 
by pattern of zeros in its derived polynomial P({2|"^}): 



I ^1 J 



p(.r 

~(^^'*)-^f))^»^p(.j('^+''), 

where Dab satisfy 



=2(0+!,) 



(4) 



Dab = Dba e Z, Daa = even. Dab > 0. (5) 

{Dab] is the set of data that specifies the local condition 
that $({zi})'s must satisfy. Such a set of data is called 
a pattern of zeros. 

We note that there are many different ways to fuse a Zi- 
variables into a z^"^ variable. The different ways of fusion 
may lead to different derived polynomials which are lin- 
early independent. Here we will impose a unique- fusion 
condition on the symmetric polynomial $({zi}): The de- 
rived •polynomials obtained from different ways of fusions 
are always linearly dependent, ie the derived polynomi- 
als form a one-dimensional linear space. In this paper, 
we will study symmetric polynomials ^{{zi}) that satisfy 
this unique- fusion condition and are characterized by the 
data Dab- 

Not all possible choices of {Dab} are consistent. Only 
certain choices of {Dab} correspond to symmetric poly- 
nomials ^{{zi}). The key is to find those {-Dafcj's that 
can be realized by some polynomials $({zi}). 

To get a feeling what a consistent set of {Dab} may look 
like, let us consider the following symmetric polynomials 
(the Laughlin state^) 



'^l/,{{z^}) = - Zj)-^ 



(6) 



where q is an even integer. Such a symmetric polynomial 
leads to the following derived polynomial 

Pi/,{{4''^}) (7) 

= { n [ 11(4"^ - } { n [ IK^^"' - ^i^'r'] } 



So the symmetric polynomial $i/q is specified by the 
pattern of zeros: 

Dab = qab, a, 5 e {1, 2, 3, • • • } 

where is a positive even integer. 

B. Sa characterization of polynomials 

There is another way to implement local conditions 
on a translation invariant symmetric polynomial ^{{zi})- 
We introduce a sequence of integers Sa, a — 0, 1, 2, • • • 
and require that the minimal total powers of zi , • • • ,Za 
in ^{{zi}) is given by ^o-^^ (Here Sq is defined to be 
0.) Thus, in addition to {Dab}, we can also use {S„} to 
characterize a symmetric polynomial. For a translation 
invariant symmetric polynomial, $(0,^2,- •• ,Zn) 7^ 0. 
Thus = 0. 

The two characterizations, {Dab} and {Sa}, are closely 
related. One way to see the relation is to put the sym- 
metric polynomial $(21 , • • • , zat) on a sphere as discussed 
in Appendix A. Let AT^ be the maximum power of zi 
in $(2:1, ■ • • , zn)- Then ^{zi. ■ • ■ , zn) can be put on a 
sphere with flux quanta and each variable Zi carries 
an angular momentum J = 

From the discussion near the end of Appendix A, 
We find that each type-a particle described by z'f'' in 

P{{z[°^^}) carries a definite angular momentum, which is 
denoted as ,/„. Since the lowest total power of zi, ■ • • , Zq 
is Sa, the minimal total quantum number for those 
variables is —aJ + Sa- Therefore the angular momentum 
of the z^°'^ variable is 

Ja = aJ-Sa- (8) 

Since z^^^ = Zi, we find that Ji = J. 

Again, according to the discussion near the end of Ap- 
pendix A, if we fuse two variables z^"'^ and z^''^ into ^;("+'') , 
the type- (a + 6) particle described by will carry an 

angular momentum 

Ja+b = Ja + Jb — Dab- (9) 

We see that Dab can be expressed in terms of Sa- 

Dab = Sa+b ~ Sa — Sb- (10) 

The conditions on Dab, (5), can be translated into the 
conditions on Sa- 



S2a = even, Sa+b > Sa + Sb- 



(11) 



From the recursive relation Ja+i ~ -Ja + Ji ^ Da.i, we 
find Sa+i = Sa + Da,i- Using = 0, we see that Sa can 
also be calculated from Dab- 



a-l 

Sa = yDbA- 



6=1 



(12) 



a<b i,j 



a i<j 



Due to the one-to-one correspondence between {Dab} 
and {Sa}, we will also call the sequence {Sa} a pattern 
of zeros. 
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C. Boson occupation characterization 

The symmetric polynomial ^{zi, - ■ ■ ,Zn) can be writ- 
ten as a sum of polynomials described by boson occupa- 
tions 

where ^{a,} is a boson occupation state with ni bosons 
occupying the orbital. Mathematically, is 
given by 

N 

^{n,}{zi,--- ,ZN) = Y.]14{iy (13) 

P i=l 

where P is a one-to-one mapping from {!,••• ,A^} — > 
{1, • • • , A''}, is the sum over all those one-to-one map- 
ping, and Zi, i = 1, 2, • • • , is a sequence of ordered integers 
such that the number of I valued liS is n/. 

What kinds of boson occupations {n;} appear in the 
above sum? Let us set z\ = in <^{{zi}). Since 
$(0,22,- •• ,2jv) 7^ due to the translation invariance, 
there must be a boson occupation {hi} in the above 
sum that contains one boson occupying the 2'"° orbital. 
Now let us assume that a boson occupies 2;'^°, and bring 
the second particle Z2 to 0, the minimal power of Z2 in 
^(O, 22,-- - ,Zn) is Dii: 

$(0, 02, • • • , Zn) ~ -22''" ^2(^23, ^4, • • • ) + o{z^''). 

So among those {hi} which have one boson occupying the 
0'^° orbital, there must an {hi} that contains a second 
boson occupying the 2'^ orbital where I2 = Dii = S2 — S1. 
Next, let us assume that two bosons occupy z'^" and 
z'^ orbitals and we bring the third particle Z3 to 0, the 
minimal power of Z3 is D21: 

P2(^3, Z4,---)^ Z^"P3{Z4, 05, • • • ) + o{z^''). 

Thus among those {hi} which have two bosons occupying 
the z^^^ and 0'^ orbitals, there must be an {hi} that 
contains a third boson occupying the z''' orbital where 
^3 = -D21 = S3 — 82- This way we can show that there 
must an {hi} such that the a**^ boson occupies the orbital 
with la = Sa — 5*0-1- Here a = 1, 2, • • • and li = 0. 
Let ni be the numbers of la = Sa — Sa-i that satisfy la = 
I. We see that the boson occupation state ^{ni}{{zi}) 
happens to be the state with its a**^ boson occupying the 
orbital z'". This allows us to show that ^{zi, ■ ■ ■ ,Zn) 
has a form 

^{Zi}) = ^in,}{{^^}) + E Cin,}^in,}{{z^}), (M) 
{ni} 

Or in other words 

(*{n,}|$)^0. (15) 



The two sequences, {Sa} and {n;}, have a one-to-one 
correspondence. We will call {ni} the boson occupation 
description of the pattern of zeros {Sa}- 

The boson occupation distributions {hi} that appear 
in the sum in (14) satisfy certain conditions. First the 
boson occupation {hi} can be described by a pattern of 
zeros {Sa}- Then the conditions on {hi} can be stated 
as Sa > Sa- So the minimal total power of zi, • • • , Za in 
^{ni}{{zi}) is Sa which is equal or bigger than Sa- 

Haldane^^ has conjectured that hi's in the expression 
(14) can be obtained from ni by one or many squeezing 
operations. A squeezing operation is a two-particle oper- 
ation that moves one particle from the orbital z'^ to z'l 
and the other from z'^ to z^'^, where h < I'l < I2 < h, 
and li + I2 = I'l + I2- We can show that, if hi is obtained 
from ni by squeezing operations, then the minimal total 
power of Zi, ■ • • , Za in ({zj}) is equal or bigger than 
Sa- This is consistent with the above discussion. 

Let Po,Ja be a projection operator acting on the state 
$ on sphere. -Po,j„ projects into the subspace where 
a particles in $ have a total angular momentum equal 
to Ja or less. We see that for a symmetric polynomial 
$(zi, ■ • • ,zn) described by a pattern of zero Sa, it satis- 
fies 

Pn,Jn ■ ■ ■ P3,J3P2,J2^{zi, ■ ■ ■ ,zn) = ^'(^i,-- • ,Zn), 
where Ja = aJ — Sa- This allows us to obtain 

PN,JM---P3,J3P2,J2^{n,}{zi,-- - ,Zn)j^O, (16) 

where n; is the boson occupation description of Sa- 

IV. THE CONSISTENT CONDITIONS ON THE 
PATTERN OF ZEROS 

For a translation invariant symmetric polynomial 
i>({zj}), the corresponding pattern of zeros {Dab} and 
{Sa} satisfies some special properties. Here we would 
like to find those properties as much as possible. Those 
properties will be called consistent conditions on the pat- 
tern of zeros. If we find all the consistent conditions, a set 
of integers {Dab} or {Sa} that satisfies those consistent 
conditions will correspond a translation invariant sym- 
metric polynomial ^{{zi}). We have already find some 
consistent conditions (11) and (16) on {Sa}- Here we 
would like to find more conditions. 



A. Concave condition 

If we fix all variables z|"'' except z["'\ then the de- 
rived polynomial Pa({z|"''}) gives us a complex function 
f{z['^^). The complex function f{z["'^) has isolated zeros 
at zf'^'s and possibly also at some other points. 

Let us move around two points z[^^ and z^^-* . The 
phase of the complex function f{z[°'^) will change by 



FIG. 1: Wa,bc obtained by moving z[ along a large loop 
around z^''' and counts the total numbers of zeros of 
/(a'"-*) in the loop. The crosses mark the zeros of f{z['^^) 
not at z['''' and z['^\ 

^'^Wa,bc where Wa,bc is an integer (see Fig. 1). Since 
fizi^^) has an order Dab zero at z^'^ and an order Dac 
zero at z[''\ the integer Wa,bc satisfy 

Wa,bc >Dab + D 

ac 

because /(^J"^) has no poles. Now let z^^ —>■ z[''^ to fuse 
into In this hmit Wa^bc becomes the order of zeros 

between 2;^"^ and zC'+'^'l: Wa,bc = -Da,6+c- Thus we obtain 
the following conditions on Dab 

Da.b+c > Dab + Dac (17) 

The concave condition (17) is equivalent to a condition 
on 5a: 

Sa+b+c + Sa + Sb + Sc > Sa+b + Sb+c + Sa+c- (18) 

We note that the Laughlin state ^i/g{{zi}) = fl^^jizi — 
ZjY saturates the above conditions: Da^b+c = Dab + Dac 

or Sa+b+c + Sa + Sb + Sc = Sa+b + Sb+c + Sa+c- 

B. Symmetry condition 

If wc fix all variables zf"'' except z['^\ Z2^\ and z^''\ 
then the derived polynomial Pa{{z'l"''^ }) gives us a com- 
plex function f{z["'\z2^\z!^'^^). Let us assume z["'\ Z2\ 
and z'^^ are very close to each other and far away from 
all other variables. f{z^\ z'f^ ,zf^) has a D*^ order zero 
as z^^ —>■ Zj^' and a D*^ order zero as zj^"-* 2:3'^'' . Thus 
Da,b+c — Dab — Dac is the number of zeros of f{z['^^) in 

the same neighborhood that are not at Zj*"^ and z^"^^ (see 
Fig. 2). 

Here we would like to assume that $({zi}) satisfies the 
unique-fusion condition. In this case, the type-o variables 

in the derived polynomials have "no shapes" and can be 
treated as points. Since other variables are far away, the 
zeros of /(zj"\ 22**^ Zg'^^) must satisfy certain symmetry 
conditions (see Fig. 2). 

We see that when a = b = c, /(zj;''\ Zj""*? ^3"'') can 
be zero only when z'f'^ zj"-* or when z'f\ z^\ and 



(a) (b) (c) 

FIG. 2: Pattern of zeros of f{z^^\z'^\z^^) (when viewed as 
a function of Zj"'): (a) Da,b+c — Db — Dc = 1, (b) Da,b+c - 
Db — Dc = 2, and (c) a = 6 = c and Da,a+a — Da — Da = 2. 
The dash lines form two equilateral triangles. The zeros that 
are not located at any variable are marked by crosses. 

Zg"'' form an equilateral triangle. Thus the zeros of 

f{z'f\z^\ Zg"'') (when viewed as a function of zj"'') that 
are marked by the crosses must appear in pairs. Wc find 
that Da^a+a — Da — Da must be even, or equivalently 

Sia -Sa = even. (19) 

C. n-cluster condition 

The structure of symmetric polynomials with infinite 

variables is very complicated and hard to manage. Here 
we would like to introduce an n-cluster condition that 
makes a polynomial with infinite variables behave more 
like a polynomial with a finite number of variables. A 
symmetric polynomial satisfies the n-cluster condition if, 
after we fuse the variables of $({zj}) into n- variable clus- 
ters, the derived polynomial 

Pe(^w,.-.,z(,"))~n(4"^-#^)' (20) 

i<j 

has a simple Laughlin form where / is a positive integer. 

To see the structure of cluster form more clearly, let 
us assume that the polynomial <I>(zi, • • • , zpf) describes a 
FQH state with filling fraction i/. This means that as a 
homogenous polynomial, the total order of the Zj, ^jv, in 
$(zi, • • • , zjv) satisfies 

SN = ^N^ + OiN). 

This motivates us to write $(zi, • • • , z^r) as 

$({Z,})=G({Z,;})$,.({Z,}) 

M{zi})=U('i-zjr\ (21) 

Here G({zi}) satisfies 

G(Azi,--- ,Azjv) = A^"G(zi,--- ,zjv), sn = 0{N). 

Note that G{{zi}) is in general not a single valued fimc- 
tion, since ^^{{zi}) is in general not single valued. But 



7 



the product of G{{zi}) and ^v{{zi}) is a single-valued 
symmetric polynomial. 

We can fuse the variables in G{{zi}) to obtain a de- 
rived function Gdz^-"^}) [just as how we obtain the de- 
rived polynomial P{{z'^°''^}) from the original symmetric 
polynomial <E>({zi})]. Similarly, we can also fuse the vari- 
ables in ^^{{zi}) to obtain a derived function $,^({2;j"^}): 

Mi^"^}) = n (^1"^ - 'fr'^' n (^1"^ - 'i^y^''- 



We also find that 



i,j;a<b 



z<j;a 



(22) 

Thus the derived polynomial P{{zl"'^}) can be expressed 
as 



P{{zl"^}) = Gi{zr})M{zr})- (23) 

(23) can be viewed as a definition of G({2;|"''}). 

Assuming ^{{zi}) has an n-cluster form, then if we 
fuse the variables of G{{zi}) into n-variable clusters, the 
derived function 

G({.1")}) = 1. 

Here we will require G{{z^'^^}) to satisfy more strict con- 
ditions 



G({z,|"^}) = G({.™}), 



(a%n)-| 



G(. 



,) = G(.. •,•••) 



(24) 



The second condition states that G(- • • ,zl"'\---) does 

not depend on zl"-* if a%n = 0. If G satisfies (24), we will 
say the corresponding symmetric polynomial $({zj}) = 
G{{zi})^v{{zi}) to have an n-cluster form. 

For a symmetric polynomial $({zi}) of an n-cluster 
form, its pattern of zeros Dab can be written as 



Dab = V ^ab + dab 



(25) 



where dab satisfy 



dab = if b%n = 

da,h-\-n — dab 



(26) 



The pattern of zeros Dab that satisfies the above condi- 
tions is said to have an n-cluster form. Note that v^^ab 
in (25) describe the pattern of zeros in the derived func- 
tion <l'i/({2j^"'}) (sec (22)) and dab describe the pattern 

of zeros in the derived function G{{z'f^}). 

Setting (a, 6) = (n, n) and {a,b) = (l,n) in (25), we 
find that 



or 



V ^n^ = even 



m 



z/ n = integer . 



mn = even. 



(27) 



Da,b+n = Da,b + am 



Let 



Sa = Sa- -^a{a - 1) 

We find that (see (10)) 

dab = Sa+b - Sa - Sb 

The cluster conditions (26) become 

Sa+n Sa = 0, 

Sa+n Sa — Sb-\-n ^6* 



(28) 



(29) 



(30) 



(31) 



Since = si = (see (8)), we find that s„+i = s„ and 

Sa+n - Sa = Sn- ThuS 

Sa+kn = kSn + Sa, a = 1, 2, • • • , 00. (32) 

This allows us to obtain Sa for any a > from 
si, S2, • • • , s„. Similarly, all the Sa's can be determined 
from Si, S2, ■■■ , Sn- 

Tfl 

Sa+kn = Sa+kn + ^(" + + ku - 1) 

= kSn + Sa + ■^{a + kn){a + kn-l) 

= Sa- ^a{a - 1) + k[Sn - ?(n - 1)] 
2n 2 

H (a + kn){a + kn — 1) 

2n 

_ , _ kik — l)nm , , . 

= Sa + kSn + — ^ + kma (33) 



The above result is actually valid for any positive integer 
a. 

It is convenient to introduce 

ha - Sa ^Sn-Sa ^ + ^ . (34) 

Prom (32), we can show that is periodic: 

use _ ISC 

"a — "-a+n- 

Since Si = 0, we see that hi'' = —Sn/n and = h^a ~ 
aifi- Prom (29), we see that Sa can be calculated from 



Sa = K - ah!f 



a{a — l)m 
2n ■ 



(35) 



(34) and (35) imply that the two sequences of numbers, 
{Sa} and {h^a}^ have a one-to-one correspondence and 
can faithfully represent each other. In this paper, we 
will use both sequences to characterize the symmetric 
polynomials. The h^a characterization turns out to have 
a close relation to the CPT description of the FQH states 
(see Section B). "-40-42 
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If Sa has the n-chistcr form (33), then the corre- 
sponding boson occupation numbers ri; have some nice 
properties. From (33), we see that Sa+n — Sa-i+n = 
Sa — Sa-i+m. Thus li = Si — Si-i satisiy k+n = k+m. 
This means that the boson occupation in the orbitals z'' 
has a periodic structure: every time we skip n boson, 
we skip m orbitals. Or in other words, if we know the 
occupation distribution of the first n bosons, the occu- 
pation distribution of second n bosons can be obtained 
from that of first n bosons by shifting the orbital index I 
by m. Thus the occupation numbers ni satisfy ni = ni%rn 
(see (52)). Also, each m orbitals contain n bosons. Due 
to the one-to-one correspondence between Sa and ni, we 



can also use no , • 
in mzi}). 



-1 to describe the pattern of zeros 



D. Translation invariance 

To study the translation invariance of the symmetric 
polynomial ^{zi, - ■ ■ , 2:jv), let us put the polynomial on a 
sphere (see Section A) and study its rotation invariance. 
In fact, in this paper, when we mention translation in- 
variance we actually mean rotation invariance on sphere. 

Let be the number of flux quanta going through the 
sphere. Then, each variable Zi in $(21, • • • , zn) carries an 
angular momentum J = What is the total angular 

momentum of ^{zi, • ■ • , zw)? In general, ^{zi, • • • , zjv) 
does not carry a definite angular momentum. So here 
we will calculate the maximum angular momentum of 
^{zi, ■ ■ ■ , Zn) from the pattern of zeros Dab- 

The maximum angular momentum is nothing but the 
angular momentum of 2:^^^ - the particle obtained by fus- 
ing all the N electrons together. The angular momentum 
of z^^) is given by (see (8) and (12)) 



Jn = J\ 



tot 



NJ 



E 

a=l 



Da,l =NJ-Sn 



(36) 



If Jm = for a symmetric polynomial ^({zj}), then 
(^{{zi}) is invariant under the 0(3) rotation of the sphere. 
In other words, <I>({zi}) is translation invariant. 

However, for an arbitrary choice of N and J, J/v is not 
zero in general. Jjv can be zero only for certain combina- 
tions of {N, J). For the filling fraction v =l/q Laughlin 

state, Jn=NJ- We find Jjv = if 



2J = N4, = qN - q. 



(37) 



This is the relation between the number of magnetic 
flux quanta, A^^, and the number of electrons. A, of the 
v = 1/q Laughlin state if the Laughlin state is to flU the 
sphere completely (which gives rise to a rotation invari- 
ant state). 

Assume that the symmetric polynomial ^{zi, • • • , Zn) 
has an n-cluster form described by the data 

(m: Sa\a=i.--- .n)- If we put the polynomial on a 
sphere, the maximum total angular momentum of 



<[>(zi,-- - ,Z]^) is given by (36). If A' = nNc, we find 
from (33) that 



SnNc — Sn+(^Nc-l)n — NcSn + 



mnNc{Nc - 1) 



Jt, 



JnNr - NrSr, 



mnNc{Nc - 1) 



(38) 



When A' = nN^ $(zi,-- - ,zn) can give rise to a 
Laughlin wave function (20) after fusing Zj's into A^^ 

2;|"^'s (see (20)). Since it is always possible to fill the 
sphere with the Laughlin state, this implies that there 
exist an integer 2 J to make Jtot = 0. Such an integer is 
given by 

'^SnNc '^Sn 



2J = Ns = 



nNr. 



+ m{Nc - 1) (39) 



This requires that 



2Sn = mod n. 



(40) 



To summarize, the n-cluster condition requires that if 
N%n = and (N0,N) satisfies (39), then the symmet- 
ric polynomial ^(zi, • • • , zat) must represent a rotation 
invariant state on sphere. The existence of such rotation 
invariant state requires Sn to satisfy (40). 



V. CONSTRUCTION OF IDEAL 
HAMILTONIANS 

We have seen that the pattern of zeros in an electron 
wave function $({2j}) can be described by a set of inte- 
gers 52, /Ss, • • • . In this section, we are going to construct 
an ideal Hamiltonian on sphere to realize such a kind of 
electron wave function as a ground state of the Hamilto- 
nian. 

On sphere, the set of integers Sa also has a very phys- 
ical meaning. For an electron system on a sphere with 
A^ flux quanta, each electron carries an orbital angu- 
lar momentum J = if the electrons are in the flrst 
Landau level. For a cluster of a electrons, the maximum 
allowed angular momentum is aJ. However, for the wave 
function ^{{zi}) described by Sa, the maximum allowed 
angular momentum is only Ja = aJ — Sa- The pat- 
tern of zeros forbids the appearance of angular momenta 
aJ — 5*0 + 1, aJ—Sa+2, - ■ - , aJ for any a-electron clusters 
in $({z.}). 

Such a condition can be easily enforced by a Hamil- 
tonian. Let Pg"'' be a projection operator that acts 

on a-electron Hilbert space. Pg^^ projects onto the 
subspace of a electrons with total angular momenta 
a J — 1, • • • , a J. Now consider the Hamiltonian^^'*^'^^ 



H. 



{Sa} 



E E 



p(«) 



(41) 



a a-electron clusters 



electron clusters 



sum over all a-electron clusters. 



where J2a- 

The wave function ^{{zi}) with a pattern of zeros de- 
scribed by Sa, if exist, will be the zero-energy ground 
state of the above Hamiltonian. 



9 



Wc note that the Hamiltonian H^s^} ^^^^ defined for 
any choice of Sa- However, for a generic choice of Sa, the 
zero-energy ground state of -ff{s„} may not be the one 
with a pattern of zeros described by Sa- This is because 
when we say the wave function ^{{zi}) has a pattern of 
zeros described by Sa, we mean two things: 

(a) The angular momenta aJ — S'a + 1, •■• a J do not ap- 
pear for any o-electron clusters in ^{{zi}); 

(b) The angular momenta aJ — Sa does appear for a- 
electron clusters in $({2j}). 

The zero-energy ground state of H{Sa} satisfies the con- 
dition (a). But sometimes, we may find that the con- 
dition (a) also implies that o J — Sa does not appear for 
a-electron clusters in ^{{zi}) for certain values of a. This 
means that the zero-energy ground state of i?{Sa} s-c- 
tually described by a pattern of zeros Sa which satisfy 
Sa > Sa- However, for a certain special set of {Sa} that 
describe the pattern of zeros of an existing symmetric 
polynomial, wc have Sa = Sa- For those Sa, the zero- 
energy ground state of -ff{s„} is described by the pattern 
of zeros of {Sa} itself. 

We have seen that, for a FQH state described by a 
pattern of zeros {S'a}, a state of a-electron clusters has a 
non-zero projection into the space 'Ha,Sai where Ti-a^Sa is 
a space with a total angular mentum aJ — Sa- However, 
different positions of other electrons may lead to differ- 
ent images in the space 'Ha,Sa- Let Ha be the subspace 
of 'Ha,s„ that is spanned by those images. In general 
Ha.Sa 7^ ^a- So, in general, the zero-energy ground state 
of the ideal Hamiltonian -ff{s„} may not be unique. In an 
attempt to construct an ideal Hamiltonian for which the 
FQH state $ is the unique ground state, we can add ad- 
ditional projection operators and introduce a new ideal 
Hamiltonian 

^{S.}=E E (Ps^J+PuJ^ (42) 

d a-electron clusters 

where is a projection operator into the space Ha 
and Ha is a subspace of Ha,Sa formed by vectors that is 
perpendicular to Ha- 

VI. SUMMARY OF GENERAL RESULTS 

In last a few sections, we have considered a subclass 
of symmetric polynomials $({2;j}) (of infinity variables) 
that satisfy (a) a unique-fusion condition (see discussion 
in Section III A), (b) an n-cluster condition (see discus- 
sion in Section IV C), and (c) translation invariance 

'^i{zi}) = H{zi-z})- 

The unique-fusion condition requires that when we fuse 
the variables Zi together to obtain new polynomials, we 
will always get the same polynomial no matter how we 
fuse the variables together. The n-cluster condition re- 
quires that if we fuse all the variables Zi into clusters of 



n variables each, the resulting polynomial of the clusters 

has the Jastrow form ni<j(-2^1"^ ~ ^j"^)'. 

We find that each translation invariant symmetric 
polynomial $({zi}) of the n-clustcr form and satisfy- 
ing the unique-fusion condition is characterized by a 
set of non negative integers {m;S2,--- ,Sn)- However, 
not all sets of non negative integers {m:S2,--- ,Sn) 
can be realized by such symmetric polynomials. The 
{m;S2,--- , Sn) that correspond to existing translation 
invariant symmetric polynomials (that satisfy the n- 
cluster and the unique-fusion conditions) must satisfy 
certain conditions. 

First m and Sn must satisfy (see (27) and (40)) 

TO > 0, mn = even 

2S„ = mod n (43) 

From TO, S'2, • • • ,Sn and Si = 0, we can determine Sa for 
any a > 1: (see (33)) 

Sa+kn = Sa + fc5„ + ^ + kma (44) 

Those Sa must satisfy (see (11) and (18)) 
A2(a, a) = even, 

A2 (o, 6) > 0, A3 (a, 6, c) > 0. (45) 

where 

A2{a,b) = Sa+b - Sa - Sb, (46) 

A3 (a, b, c) = Sa+b+c - Sa+b - Sb+c - Sa+c + Sa + Sb + Sc- 

The Sa's also satisfy another condition which is harder 

to describe. To describe the new condition, we first note 
that the sequence {Sa} can be encoded by another se- 
quence of non-negative integers ni, I = 0,1, ■■■ . To 
obtain n; from Sa, wc introduce la = Sa — Sa-i for 
a = 1,2,- ■•. Then n; is the number of /a's that sat- 
isfy la = I- The two sequences, {Sa} and {n;}, have a 
one-to-one correspondence and can faithfully represent 
each other. The number m can be regarded as the boson 
occupation number that was used to characterize FQH 
states in the thin cylinder limit. ^^"^^ n; is also used to 
label Jack polynomials that describe FQH states. ^^'^^ 

Now let us introduce 2J + 1 orbitals \m,z), = 
—J,—J + l,--- , J — 1, J, which form a representation 
of SU{2) with an angular momentum J. (Here 2 J is an 
integer.) We can create a many-boson state |{n;}) by 
putting ni bosons into the ruz = I — J orbitals. Then Sa 
must be such that (see (16)) 

Pn.nj-s^ ■ ■ ■ Pi,zj-s,P2.,2j-s,\{ni}) ^ 0, (47) 

where Pa,J„ is a projection operator that projects into 
the subspace where any a particles have a total angular 
momentum equal to Ja or less and N is the number of 
particles in |{n(}). 
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We will call (m; 52, • • • , Sn) an S-vector and denote it 

as 

S = (m;S'2,-- - ,Sn) 

We find that translation invariant symmetric polynomials 
(that satisfy the n-cluster and the unique-fusion condi- 
tions) are labeled by the S- vectors that satisfy (43), (45), 
and (47). 

We would like to stress that (43), (45), and (47) are 
only necessary conditions for (m; 5*2, • • • , Sn) to describe 
a translation invariant symmetric polynomial of the n- 
cluster form and satisfying the unique fusion condition. 
Wc do not know if those conditions arc sufficient or not. 
Some 5- vectors that satisfy (43), (45), and (47) may not 
correspond to an existing symmetric polynomial. Also, 
there may be more than one symmetric polynomials that 
are described by the same 5- vector that satisfy (43), (45), 
and (47). Each such symmetric polynomial corresponds 
to a FQH state. By solving (43), (45), and (47), we can 
obtain (m; ^2, • • • , 5'n)'s that correspond to the Laughlin 
states, the PfFafian state, the parafermion states,"^® and 
many new non-Abelian states. 

We also obtained some additional results. Our numer- 
ical studies of the equations (43) and (45) suggest that 
all solutions of the equations satisfy 



where 



hT = hi 



Lsc _ c ^^"^ 
"a — Sa I — , 

n 2 2n 



(48) 



(49) 



out to be the scaling dimensions of those simple-current 
operators. Since A^(a,b,c) only depends on h^^ 

A3(a,6,c) = /I- ^-/i- ^ + + + 

A3 (a, 6, c) > is a property of the simple-current CFT. 

The condition (47) is hard to check. So let us just 
consider (43) and (45). One class of solutions of (43) and 
(45) is given by 

j^sc ^ ^z^ = {a%n)[n-{a%n)] ^ 



This class of solutions corresponds to the Z„ parafermion 
CFT which is generated by simple-current operators t/'a 
that satisfy an Abclian fusion rule i/'aV'b — i'a+b and 
V'o = V'n = 1- The scaling dimensions of ipa is given 
by the above /if" . The parafermion states introduced in 
Ref. 19 are related to such a class of solutions. 

A more general class of solutions of (43) and (45) corre- 
sponds to generalized parafermion CFTs. A generalized 
parafermion CFT is generated simple-current operators 
that have the following dimensions 

sc _ f^z^ _ {ka%n)[n - {ka%n)] 
n 

Those solutions represent a new class of non-Abelian 
FQH states, which will be called generalized parafermion 
states. 

It turns out that all the solutions of (43) and (45) are 
closely related to parafermion CFTs, ie a solution 
satisfies 



although we cannot derive the equation (48) analytically. 
Such a relation implies that 



ki 



(50) 



n — 2a 

Sn—a = 5a H Sn- 



n 



h^^s also satisfy 



tSC _ LSC LSi 



> 0. 



where a%n = a mod n. 

From (49), wc find that (hf, ■■■ , /i^f ) and (^2, • • • , 5„) 
have a one-to-one correspondence. They can faithfully 
represent each other. Due to the one-to-one relation be- 
tween Sa and h^^, we can also use n, m, and hf^, - ■ ■ , h'^ 
to characterize the pattern of zero in the symmetric poly- 
nomial ^{{zi}). We will package the data in the form 



h={-;hr 



and call h an /i-vector. We see that patterns of zeros 
in a symmetric polynomial can also be described by the 
/i- vectors. 

Each symmetric polynomial described by the pattern 
of zeros {5a} is related to a conformal field theory (CFT) 

generated by simple-current operators which have an 
Abelian fusion rule (see section B). {hf', ■ ■ ■ ,h^n) turn 



where fcj's are positive or negative integers and /Iq ''s 
are the scaling dimensions of the parafermion operators 
in some parafermion CFT's labeled by i. They are given 
by 



h 



PF 



/la" 



for certain integers k and n', where n' is a factor of n. 

If h^^ = h^^ , then the solution corresponds to an ex- 
isting symmetric polynomial generated by a (generalized) 
parafermion CFT. If h^^ = \h^^ , then the solution cor- 
responds to the square root of a symmetric polynomial 
generated by a (generalized) parafermion CFT. There- 
fore, the later solution does not correspond to any exist- 
ing symmetric polynomials. Numerical experiments sug- 
gest that the later cases always have A3(a, 6, c) = odd 
for some a, b, and c. This motivates us to introduce the 
following new condition 



A3 (a, 6, c) = even 



(51) 



to exclude those illegal cases. The conditions (43), (45), 
and (51) provide an easy way to obtain 5a 's that may 
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correspond to existing symmetric polynomials. The new 
condition (51) is a generalization of necessary conditions 
A3(a, a,a) = even (see (19)), (B8), and Sa 7^ 1. 

Our numerical studies suggest that the solutions of the 
equations (43), (45) and (51) give rise to h^^ that satisfy 



B. Products of symmetric polynomials 



hf = ^hhl^^ mod 2 



Those solutions also have the properties that m = even 
and Sa = even. 

We also find that, for Sa satisfying (43) and (45), the 
corresponding n; is a periodic function of I (for / > 0) 
with a period m: 



The n;'s satisfy 

2J 



2.J 



^ n; = nNc, ^(Z - J)ni = 0, 



(52) 



(53) 



1=0 1=0 
for any (J, Nc) satisfying 

2 J = — + m(N^ - 1) 
n 



(54) 



VII. GENERAL STRUCTURE OF THE 
SOLUTIONS 

The S'-vectors that satisfy (43) and (45) has some gen- 
eral properties. In this section we will discuss those prop- 
erties. 



A. n-cluster polynomial as /tn-cluster polynomial 

Let P{{z^'^^}) be a derived symmetric polynomial of 
n-cluster form described by {m; S2, Sn). Prom (26), 

we see that P({4'^^}) can also be viewed as a symmetric 
polynomial of Kn-cluster form where k is a positive inte- 
ger. When viewed as a Kn-cluster polynomial, P{{zl"'^}) 
is described by {kiti, S2, Sku), where 5„+i,-- - 
are obtained from (m;S'2, ■■■,Sn) through (33). 

The the filling fraction u = l/q Laughlin state $1/^ = 
ni<j('^i — ZjY has a 1-cluster form. Thus $1/^ can also 
be viewed as an n-cluster polynomial for any positive n. 
When viewed as an n-cluster polynomial, the u = 1/q 
Laughlin state is described by 

(rn; A2, • • ■ , ^n) = [nq; • • • , )■ 

Such a v = 1/q Laughlin state always appears as a solu- 
tion of (43) and (45) for any n. 



Let P({z-°^}) and P'{{zl°''}) be two derived symmetric 
polynomials of n-cluster form described by (m; 5*2, Sn) 
and (m'; 5*2, 5*^) respectively. Then their product 

P({4"^}) = P({4"^})^"({4"^}) is also a symmetric 
polynomial of n-cluster form. P{{zl"'^}) is described by 

(m, ^2, • • • , Sn) = {m + m'; S2+S2,--- ,Sn + S'J. 

This is because the pattern of zeros of Pdzl"-*) is re- 
lated to the patterns of zeros of P{{z^°'^}) and P'{{z'i°'^}) 
through 



D. 



ah 



D. 



ah 



D' 



ah 



Also the relation between (m; 6*2, Sn) and Dab is lin- 
ear (see (10), (12) and (33)). Therefore, if two S^-vectors, 
S and S' , describe two existing symmetric polynomials, 
then their sum S = S + S' also describes an existing 
symmetric polynomial, whose fillings fractions are recip- 
rocally additive. 

Indeed, the solutions of (43) and (45) have a struc- 
ture that is consistent with the above result. We note 
that (43) and (45) arc linear in the S'-vector S = 
{m;S2,-' ,Sn)- Thus if Si and ^2 are two solutions 
of (43) and (45), then 



S = kiSi+k2S2 



(55) 



is also a solution for any non- negative integers ki and ^2. 
Therefore we can divide the solutions of (43) and (45) 
into two classes: primitive solutions and non-primitive 
solutions. The primitive solutions are those that cannot 
be written as a sum of two other solutions. All solutions 
of (43) and (45) arc linear combinations of primitive so- 
lutions with non-negative integral coefficients. 

As an application of the product rule, let us consider 
a symmetric polynomial of n-cluster form <I>({zj}) which 
is described by (m; ^2, • • • ,Sn)- We can construct a new 
symmetric polynomial of n-cluster form from $({2;,}): 



H{z.}) = H{z.})l[{ 

i<j 



where q is even. The symmetric polynomial $({2;^}) is 
described by 

/ ~ ;^ ;^ \ / r, qnin — 1) , 
(m; ^2, • • • , Sn) = [m + nq; S2 + q, ■ ■ ■ , Sn -\ ^ ) 



VIII. SOME EXAMPLES 

In this section, we will give some examples of 
symmetric polynomials described by the S'-vector 
(m; 52, • • • , Sn) that satisfy (43), (45), and (51). 
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1 CEises 



and 



If n = 1, the different patterns of zeros are character- 
ized by an even integer m. We find Sa = ma{a — l)/2 
and Dab = mab. Each even m corresponds to a = 1 /m 
Laughlin state 



^<3 

Wc have introduced three equivalent ways to describe 
a pattern of zeros: the S-vector (m; S'2, • • • , Sn), the h- 
vector (^; /if'^, • • • , /i^f ), and the boson occupation num- 
ber n/: (no, • • • , rim-i)- For the v = 1/m Laughlin state 
those data are given by 



{no, ■ 



,m 



S = (m;), 
{-■hf)^{m-Q), 

.i) = (l,0,---,0). 



n 

1 ^m- 



B. n = 2 cEises 

If n = 2, the different patterns of zeros are charac- 
terized by two integers m,S2- The following two sets of 
m, ^2 are the primitive solutions of (43) and (45): 

(m;52) = (l;0), (m'; 5^) = (4; 2) 

Let us discuss the solution {m; S2) = (1;0) in more 
detail. The corresponding boson occupation numbers are 

(no,ni,---) = (2,2,2,. ••) 

where there are two bosons occupying each orbital. Let 
us check condition (47) for the J = 1/2 case where there 
arc only two orbitals. This leads to a state |{2,2}) with 
4 bosons described by the wave function 

'^{2,2} = Z1Z2 + 2:12:3 + ZlZi + 22Z3 + Z2Zi + Z3Z4 

On sphere the above wave function becomes (see Section 
A) 

*{2,2} = S[viV2U3U4\. 

where <S is the symmetrization operator. Since ^4 = 2, 

we find that J4 = 4J — ^4 = and Pa^j^ is a projection 
into the subspace with vanishing total angular momen- 
tum. A direct calculation reveals that ^4, 74^12 2} ~ ^• 
Thus (m; ^2) = (1; 0) does not satisfy the condition (47), 
and docs not correspond to any translation invariant 
symmetric polynomial. 

Now let consider m, ^2 that satisfy a new condition 
(51) in addition to (43), (45). The following two sets of 
m, 52 are the primitive solutions: 



l.z^: (m;52) = (2;0), 
(^;/if,-..,/i-) = (^;i 0), 
(no, • • • ,nm-i) = (2,0) 



1/2 



(m;52) = (4;2) 



771 4 

(-;/.f,-..,/i-) = (-;0,0), 



("0,- 



, Tin 



.1) = (1,0,1,0) 



Here we also listed the corresponding /i-vector h = 
{■^',hi^,- " jh^) and the boson occupation numbers 
{no, ■ ■ ■ ,nm-i)- 

Let us discuss the solution (rn; 5*2) = (2; 0) in more 
detail We find (51,52,53,54) = (0,0,2,4) and 



(Dn 


Dl2\ 






\D21 


D22 







which means that we will have no zero if we bring two 
particles together and a second order zero if we bring a 
third particle to a two-particle cluster. Such a pattern of 
zeros describes the following translation invariant sym- 
metric polynomial 



■,Z: 



(M)=^( 



1 



1 



Zl - Z2 Z3- Zi 



■)W^Zi-zj) (56) 



which is the filling- fraction 1^=1 bosonic Pfaffian state. -'^^ 
Here A is the antisymmetrization operator. The Pfaffian 
state can be written as a correlation of the following op- 
erator in a CFT 



V,{z) = i>{z)e 



l<t>(z) 



where fl}{z) is the Majorana fermion operator in the Ising 
CFT (which is also the Z2 parafermion CFT). The hf = 
1/2 in the /i- vector is the scaling dimension of ij). 

The other solution (m; 52) = (4; 2) gives rise to the 
following Dab' 



Dn D12 
D21 D22 



2 4 

4 8 



It describes the symmetric polynomial 

*i/2({^i}) = n^^' ~ 



(57) 



which is the filling-fraction i> = 1/2 bosonic Laughlin 
state. The = 1/2 Laughlin state can be written as 
a correlation of the following operator in the Gaussian 
model (or [7(1) CFT) 

Ve{z) = e'^'^^'l 

We note that the u = 1/2 bosonic Laughlin state is 

characterized by a pattern of boson occupation numbers 
(n;) = (1, 0, 1, 0, 1, 0, • • • ), and the v = 1 bosonic PfaSian 
state characterized by {ni) = (2, 0, 2, 0, 2, 0, ■ • • ). Those 
patterns match the boson occupation distributions of the 
two states in the thin cylinder limit. '^^"'^^ This appears to 
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be a general result: the ni that characterize a symmet- 
ric polynomial correspond to one of the boson occupation 
distribution of the same state in the thin cylinder limit. 
Or more precisely: the ni that characterize a symmetric 
polynomial correspond to the boson occupation distribu- 
tion of the same state in the thin sphere limit.'^^ 

The solution (m; S2) = (4; 0) = 2 x (2; 0) gives rise to 
the following Dab'- 



(58) 



It describes a symmetric polynomial which is the square 
of $i 





Dl2\ 


-I 




\D2l 


D22 







1 1 

Zl - Z2 Z3- Zi 



(59) 



Let us consider another translation invariant symmet- 
ric polynomial 



^dw{{Zi}) = S 



1 



1 



{Zi - Z2Y (z-i - Z4)^ 



i<j 

(60) 

where S is the symmetrization operator. We note that 
^dw{{zi}) and $22X2 li^s the same pattern of zeros 

given by (58). In the following, we would like to show 
that 

^dw{{Zi}) oc ^z^zA{zi})- 

We first note that ^z^z-iiizi}) can be written as a 
correlation of the following operator in a CFT 

K(^)=Ai(z)A2(z)e*^(^) 

where Ai(z) is the Majorana fermion operator in a 
Ising CFT and X2{z) is the Majorana fermion opera- 
tor in another Ising CFT. Thus the operator 14 (z) = 
Ai(2)A2(2)e*'^^^) is an operator in Ising x Ising xJ7(l) 
CFT. 

The state ^dw{{zi}) can also be written as a correla- 
tion of the following operator in a CFT 

V,{z) = d,4>{z)e'*^'^ 

where (p{z) is the field in a second U{1) CFT. Thus the 
operator Ve{z) = dz(f>{z)e^'^^^^ is an operator in U{1) x 
i7(l) CFT. 

From the bosonization of the Isingxising CFT, one can 
show that the Isingxising CFT is equivalent to the C/(l) 
CFT and A1A2 has the same A^-body correlation func- 
tions as d;,4>{z). Thus ^dw{{zi}) = ^Z2Z2{{zi})- Nu- 
merical calculations have suggested that ^dw{{zi}) has 
gapless excitation and is unstable. 



Next, let us consider the following two polynomials 
«(H.).({^.}) = *1^.({^<)) (61) 

1<J 



and 



{zi - Z2Y [zs - Zi)l 



■)W^^-z,Y, 

(62) 

where Sq = S when q = even and Sq = A when g = odd. 
The two symmetric polynomials have the same pattern 
of zeros Dab- But when q > 2, the two polynomials are 
different. Those polynomials provide us examples that 
there can be more than one polynomials that have the 
same pattern of zeros. 

C. n = 3 cases 

If n = 3, the different patterns of zeros are character- 
ized by three integers m, 82,83. The following two sets 
of m, ^2, 53 are the primitive solutions of (43) and (45): 

^s.z^: (m; 82,83) = {2; 0,0), 

(-;/.r,...,/.-) = (|;|,|,0), 
n 3 3 3 



(no,-- - ,nm-i) = (3,0) 



and 



$1/ 



2 • 



(m;52,53) = (6;2,6) 



{no,' 



,/i-) = (-; 0,0,0), 
Im-l) = (1,0,1,0,1,0) 



When n = odd, we find that the solutions of (43) and 
(45) automatically satisfy (51). 

From the /i-vector of the solution (m;52,S'3) = 
(2; 0, 0), we find that the corresponding polynomial ^3 .^^ 
describes the Z3 Rcad-Rczayi parafermion state, since 
= hf' = 2/3 in the /i- vector match the scaling 
dimensions of the simple-current operators in the Z3 
parafermion CFT. Such a state has a filling- fraction 
v = 3/2. Here we have been using $^:z„ to denote 
a Zn parafermion state. We will follow such a con- 
vention in the rest of this paper. The second solution 
(m; 52,53) = (6; 2, 6) describes the u = 1/2 Laughlin 
state. 



D. n = 4 cEises 

When n = 4, the different patterns of zeros are charac- 
terized by four integers 82, 83, 84. The primitive so- 
lutions of (43) and (45) are given by the following three 



14 



sets of m, 82,83,84: 

{m; 82, 83, 84) = {1; 0,0,0), (2; 0,1, 2), (8,2,6,12). 

The solution S = (2; 0,1,2) is the same as solution 

S = (1;0) for the n = 2 case {ie the two solutions give 
rise to the same sequence {8a}, a = 1, 2, 3, • • • ). Such a 
solution does not satisfy (47) as shown in section VIII B. 

The sohition S = (1; 0,0,0) does not satisfy (47) ei- 
ther. Let us check the condition (47) for ,/ = 1/2 case 
where there are only two orbitals. This leads to a state 
I {4, 4}) with 8 bosons. On sphere, such a state is given 
by 



{4,4} 



■ S[viV2V3V4U5UeUrU8]. 



Since 8s = 4, we find that = 8J — 8^ = and Ps.Jg 
is a projection into the subspace with vanishing total 
angular momentum. Explicit calculation shows that the 
state 4} bas a vanishing projection onto the Jtot = 
subspace. 

So we consider the solutions of (43), (45), and (51) to 
exclude those invalid cases. The primitive solutions of 
(43), (45), and (51) are 



$4 



:{m;82,--- , 5„) = (2; 0, 0, 0), 



(no,-- - ,n 
: {m;82,--- 

{no, - - 



Im-l) = (4,0), 

,5„) = (4;0,2,4), 

,K) = {l;l,0,l,0), 

im-i) = (2,0,2,0), 



(63) 



(64) 



$1/2 : (m;52,-- - , 5„) = (8; 2, 6, 12) 
(^;/if,---,/i-) = (^;0,0,0,0), 
(no,-- - ,n„_i) = (1,0,1,0,1,0,1,0). 



Among those three primitive solutions, only ^^-Zi 
new. From the /i- vector in (63), we find that the solution 
(m; ^2, 53, 84) = (2; 0, 0, 0) describes the Z4 parafermion 
state $4.7 with 1/ = 2. 

2 1-^4 

The solution (m; 82, S'3, >S'4) = (4; 0, 2, 4) is the same as 
the = 1 bosonic Pfafhan state (to;S'2) = (2;0) dis- 
cussed before. So the /i-vectors of the two solutions, 
(^;ftf,... ,/if) = (|;i,0,i,0) and {f;hf,hf) = 
(|;|,0), characterize the same state. In fact, the 
repeated (5,0) pattern in {hl'^ , ■ ■ ■ , h^) implies that 
{f;hf,---,hf) = (|;i,0, i,0) can be reduced to 
{^■,hf,hf) = (|;5,0). Also the two solutions lead to 
the same pattern of boson occupation numbers: (n;) = 
(2, 0, 2, 0, 2, 0, 2, 0, ■ • • ), again indicating that the two so- 
lutions describe the same state. 



The solution (to; S'2, S'3, S'4) = (8; 2, 6, 12) describes 
the v = 1/2 Laughlin state ^1/2 characterized by 
{^■,hf,---,hf) = (|; 0,0, 0,0). We have seen that 
hi/2 = (f;0), h2/4 = (|;0,0), /i3/6 = (f; 0,0,0), and 
— (f; 0,0,0,0) all describe the same v = 1/2 
Laughlin state $1/2- All the above solutions share the 
same pattern of boson occupation numbers: (n;) = 
(1,0,1, 0,1, 0,1, 0,---) 

The product state ^'4.2^$2.^^ described by 

*|;Z,*i;Z2 : (m;52,-- - ,S„) = (6;0,2,4) 

0'^r,-,/^?.^)-(||i,|o), 

(no,-- - ,nm-i) = (2,0,2,0,0,0). 

is a possible stable = 2/3 FQH state. Note that the 
above /i- vector is the sum of the /i- vectors of the Z2 and 
Z4 parafermion states. 



E. n = 5 cases 

When n = 5, the conditions (43) and (45) have the 
following three sets of primitive solutions: 



;4 



(to;52,- 


- - , Sn) 


= (2; 0,0, 0,0), 


—;hf,- 

n 




= (2.1 M lo)- 

^5'5'5'5'5' ^' 


{no, ■ ■ ■ 


,nm-i) 


= (5,0). 


{m;82,- 




= (8; 0,2, 6, 10), 


n 




8 6 4 4 6 

"^5'5'5'5'5' ^' 


{no, ■ ■ ■ 


,nm-i) 


= (2,0,1,0,2,0,0,0) 



$1/2 : (m;6'2, 
("0, ■ 



,5„) = (10;2,6,12,20), 

,/i^^) = (^; 0,0, 0,0,0). 
5 

Irn-l) = (1,0,1,0,1,0,1,0,1,0). 



All other solutions are linear combinations of the above 
three solutions. 

{i^;hf,---,hf) = (|;|,|,|,|,0) describes a Z5 
parafermion state ^i-z^ studied by Read and Rezayi.^^ 
(f ;/if,.-- ,/ir) ="(|;|,|,|,|,0) describes a new 
parafermion state $5 „(2) with v = 5/8. The third state 

i^;hf, • • • , hf) = (^; 0, 0, 0, 0, 0) describes the u = 1/2 
Laughlin state $1/2- 

The Z5 parafermion state ^^-z^ can be expressed as 
a correlation of simple-current operators ijji in the Z^ 
parafermion CFT. tpi has a scaling dimension of hf^ = 
4/5. The new parafermion state $ 5 „(2) can be expressed 
as a correlation of simple-current operators ^J2 in the Z^ 
parafermion CFT. ip2 has a scaling dimension of /il'^ = 
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6/5. In general, the simple-current operator (pi of a Zn 
parafermion CFT has a scaling dimension 



hf = 



l{n-l) 



Here we have been using to denote a generalized 

Zn parafermion state. We will follow such a convention 
in the rest of this paper. 



G. n = 7 cases 

When n = 7, the conditions (43) and (45) have the 
following five sets of primitive solutions: 

$i;z, : (m; 52, • • • , Sn) = (2; 0, 0, 0, 0, 0, 0), 

(H^.^sc .sc^ .2 6 10 12 12 10 6 
n 7777777 
(no,-- - ,nm-i) = (7,0); 



F. n = 6 CEises 

When n = 6, the conditions (43), (45), and (51) have 
the following four sets of primitive solutions: 



$,.^(2) : (m;S2,-- - ,5„) = (8; 0, 0, 2, 6, 10, 14), 



dH-h^o ^sc^_.8.10 12 6 6 12 10 
n 7' 7' 7'7'7' 7' 7 '' 

(no,-- - ,n™_i) = (3,0,1,0,3,0,0,0); 



(m; 52 



,5„) = (2;0,0,0,0,0), 



'''"^-^6'6'3'2'3'6' ^' 

(«0,-- - ,rim-l) = (6,0). 

$2.^2, $3 5^3, and $1/2- Three of the four primitive so- 
lutions have been discussed before and only one solution, 
$6.^^, is new. The $6.^^ state is the Zq parafermion 

state. ^^ 22 ^^"^ ^i-z the Z2 and Z^ parafermion 
states discussed before. 

Using ^2.^^, ^^"^ *^|-26' '^'^''^ construct some 

interesting and possibly stable composite states: 



^|;Z3*f;Ze ^ K ^2, • • " , 5„) = (6; 0, 0, 2, 4, 6), 

^n'^' '''"^"^6'2' '2' '2 
(no,-- - ,n„_i) = (3,0,3,0,0,0); 



:(m;52,-- - , 5„) = (8; 0, 2, 4, 8, 12), 

,/i„)-(g,3,3,2,3,3,U), 
(no,-- - ,n„_i) = (2,0,2,0,2,0,0,0); 



i>2.7„$! 



;23 



: (m;52,-- - , 5„) = (10; 0, 2, 6, 12, 18), 

fH^.h^o ... ;,sc^ .10.7 2 12 7 

'''"^-^ '6'3'2'3'6' ^' 
(no, - • - , n„_i) = (2, 0, 1, 0, 1, 0, 2, 0, 0, 0); 



*f;Z2^i;Z3^f;Z6 ■■ K^2,-- - , 5„) = (12; 0, 2, 6, 12, 18), 
(^;/if,-- - ,/i-) = (^;2,2,2,2,2,0), 
(no,--- ,n„_i) = (2,0,1,0,1,0,2,0,0,0,0,0). 

The filling fractions of those states are given by u = n/m. 



7 .7(3) 



: (m; 52, - - - , 5„) = (18; 0, 4, 10, 18, 30, 42), 



18 12 6 10 10 6 12 



(111- ^sc ... f^SC\ _ _ _ _ _ 0\ 

(no, - - - , n„_i) = (2, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0); 



$i;C, : (m; 52, - - - , 5„) = (14; 0, 2, 6, 12, 20, 28), (65) 

(^;/zr,--- ,C) = (y;2,2,2,2,2,2,0), 
(no, - - - , n„_i) = (2, 0, 1, 0, 1, 0, 1, 0, 2, 0, 0, 0, 0, 0); 

and $1/2- 

^|;Z7 is the Zr parafermion state which can be ex- 
pressed as a correlation of simple-current operators ipi 
in the Zj parafermion CFT. ipi has a scaling dimen- 
sion of /if = 6/7. $7.„(2) and $ 



7 .7(3) 



are two new 



Zj parafermion states. <I>7 (2) can bo expressed as a 

correlation of simple-current operators ip2 while $ 7 (3) 

can be expressed as a correlation of simple-current op- 
erators tps in the Z^ parafermion CFT. ip2 has a scaling 
dimension of /il'^ = 10/7 and Vs has a scaling dimension 
of hf = 12/7. 

Let us discuss the state $ 7 in more detail. $ 7 

14,07 T4''^7 

has a form 

^J,;cA{^i}) = GcA{zi})l[(^i-zjf 

where the pattern of the zeros (or poles) for Gcj{{zi}) is 
given by (see (25)) 

/-2 -2 -2 -2 -2 -4 0\ 
-2 -2 -2 -2 -4 -2 
-2 -2 -2 -4 -2 -2 

(dob) = -2 -2 -4 -2 -2 -2 

-2 -4 -2 -2 -2 -2 
-4 -2 -2 -2 -2 -2 

Vooooooo/ 

where a, 5 = 1, - • - , 7. It implies that there is a second 
order pole as an a-cluster approaches a ^-cluster if a-|-& < 
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FIG. 3: The graph that represents fcj- Each Une between 
the i**" dot and the j*^ dot represent a factor l/{zi — Zj). 



7 and a fourth order pole as an a-clustcr approaehes a b- 
cluster if a + 6 = 7. Such a pattern of poles is reproduced 

by 



where 



fCn{zi, ■ ■■ ,Zn) = 



Zi — Z2 Z2 — Z3 Zn-1 — Zn Zn — Z\ 



To confirm such a result, we note that, for a = 2, ■ • • ,7, 
the minimal total powers of a variables in /^^ are with 
(S2,---,S7) = (-2,-4,-6.-8.-10,-14Hscc (29)). 
The minimal total powers of a variables in W{zi — Zj)"^ 
are with (53, • • • , S7) = (2, 6, 12, 20, 30, 42). Thus the 
minimal total powers of a variables in ^j_.c^ are given 
by Sa = sa + §a- {S2, ■■■,87) = (0, 2, 6, 12*20, 28) which 
agrees with (65). 

The function fcj{zi, ■ ■ ■ ,27) can be represented graph- 
ically as in Fig. 3. In such a graphic representation, the 
maximum total order of the poles of a variables is the 
maximum number of lines that connect a dots. Note 
that the maximum total order of the poles of a variables 
is the negative of the minimal total power of zeros of a 
variables. 

In fact, for any n, we have a state ^j2.-c„ described 
by {hf, ■■■ , /i^'Li, h""^) = (2, • • • , 2, 0). The exphcit wave 
function is given by 



Zi Zj^ X 



^^■,c. = n( 

i<j 



> Zn)fcni^n+l, • • • , Z2n) ' 



H. n = 8 cases 

When n = 8, the conditions (43), (45), and (51) have 
the following six sets of primitive solutions: 

^|;Zs • ^2, • • • , Sn) = (2; 0, 0, 0, 0, 0, 0, 0), 

(Hl-h^c... ^sc. .2.7 3 15 15 3 7 

(no,-- - ,nm-i) = (8,0); 





(a) 



(b) 



FIG. 4: The graphs that represent /cg/Za- Ea-ch line between 
the i"^ dot and the j**^ dot represent a factor l/{zi — Zj). After 
the antisymmetrization, /cg/Za from graph (a) gives rise to 
a non-zero antisymmetric function, while /cg/Zj from graph 
(b) gives rise to vanishing antisymmetric function. 



$^.„(3) : (m; ^2, - • • , Sn) = (18; 0, 2, 8, 14, 24, 36, 48), 

18 '•^8 

(H^.h^o ... ,,sc^ .18.15 3 7 7 3 15 
[^,n, , ,A„j-i g , 8 '2'8' '8'2' 8 ' 

(no, - - - , n„_i) = (2, 0, 1, 0, 0, 0, 2, 0, 0, 0, 1, 0, 2, 0, 0, 0, 0, 0); 



■■ (m; S'2 



(8; 0,0, 2, 4, 8, 12,16), 



(^;'.r,-,'.?) = (||2,|24,2,|o), 

(no,-- - ,n„_i) = (3,0,2,0,3,0,0,0); 

*4;Z4' *i;Z2' and $1/2. 

is the parafermion state which can be ex- 
pressed as a correlation of simple-current operators V^i 
in the Zg, parafermion CFT. ipi has a scaling dimension 
of hf^ = 7/8. „(3) is a new Zg parafermion state. 

IS '^8 



can be expressed as a correlation of simple- 



current operators ■i/'s in the Zs parafermion CFT. ip^ has 
a scaling dimension of h?^ = 15/8. 

Let us discuss the state $8.^ in more detail. We 
note that the /i-vector for the ^s.-c8/Z2 state is the dif- 
ference of the h- vectors of the $ 8 state and the $2.7 
state: 

,3 „ 3 „ 3 „ 3 
(-,2,-,2,-,2,-,0) 
^2 2 2 2 ^ 

= (2, 2, 2, 2, 2, 2, 2, 0) - (i, 0, ^ 0, i, 0, i, 0). 
^|;Cs/Z2 bas a form 

*|;C8/Z2({^i}) = Gcs/Z^iiZi}) YliZi - Zj) 

where the minimal total power of a variables in 
Gc8/Z2{{zi}) is given by Sa with (s2,---,S8) = 
(-1, -3, -4, -6, -7, -9, -12). We find that 

Gcs/Z2{{Zi}) = A fcs/Z2{zi,- ■ • ,Z8)fCs/Z2{Z9,--- , zw) • 

where A is the antisymmetrization operator and the func- 
tion fcs/Z2{zi, - - - , -Zs) is represented by Fig. 4(a). 
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FIG. 5: The graph that represents fdiz-i- Each hne between 
the i**" dot and the j**^ dot represent a factor — Zj). 





FIG. 6: The graphs that represent ^c^jz-i- Each hne between 
the t'^ dot and the dot represent a factor l/(zi — Zj). After 
the antisymmetrization, fc(i/Z2 from graph (a) gives rise to 
a non-zero antisymmetric function, while fc^/z^ from graph 
(b) vanishes. 



where the function /C4/Z2 (-^ii " ' '-^4) represented by 
Fig. 5. In fact, fc^/z^{zi,- ■ ■ ,Zi) = ni<i<j<4 i"^ is 
the only function whose total order of poles for 2-, 3-, and 
4-particle clusters are given by 1, 3, and 6 respectively. 
Such a state is studied recently by Yue Yu.^^ 

Another interesting state is ^s..Cs/Z2 which has the 
same pattern of zeros as a composite state of and Zq 
parafermion states: 



*f ;C6/Z2 lias a form 



^hCe/Z^iizi}) = Gc,/zA{zi})\{{z^ - Zj) 

where the minimal total power of a variables in 
Gce/Z2{{zt}) is given by Sa with (s2,---,S6) = 
(-1, -3, -4, -6, -9). We find that 

Gca/z^iizi}) = -A fcs/z^izu- ■ ■ ,^6)/c6/Z2(^7, • • • ,2^12)- 

where the function /C6/Z2 (-^ij • • ' is represented by 
Fig. 6(a). 



I. n = 9 cases 



We would like to mention that a state simpler than 
^|;Cs/Z2 is ^|;C4/Z2 ^^^^ the same pattern of zeros 
as a composite state of Z^ parafermion state ^^-Zi'- 

4,04/^2 2''^4 2''^4 

Here « means to have the same pattern of zeros. 
*|;C4/Z2 has a form 

^|;C4/Z2({^i}) = Gc,/zA{Zi}) \{{Zi - Zj) 

where the minimal total power of a variables in 
Gc4/Z2{{zi}) is given by with (s2,---,S4) = 
(-1,-3,-6). We find that 

GC4./Z2 {{zi}) = A /C4/Z2 (^1, • • • > ^4)/c4/Z2 {z5, ■■■ , ^s) • • 



When n = 9, the conditions (43), (45), and (51) have 
the following six sets of primitive solutions: 



^%;z, ■■ (w; ^2, • • • , Sn) = (2; 0, 0, 0, 0, 0, 0, 0, 0), 
( 

(no,-- - ,n„_i) = (9,0); 



f!!^./,sc ... /,sc^ .2.8 14 2 20 20 14 8 
^n'^' ' "^~^9'9' 9' ' 9' 9' ' 9'9' ^' 



$2.^(2) : (m; ^2, - - - , 5„) = (8; 0, 0, 0, 2, 6, 10, 14, 18), 
(no,-- - ,n„_i) = (4,0,1,0,4,0,0,0); 



i;Z, 



(4) 



:{m;S2,---,Sn) = (32; 0, 6, 14, 26, 42, 60, 84, 108), 

20 

y 

(no, - • • , n„_i) = (2, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0); 



n 



32 20 8 14 14 8 
' ~ ^Y'T' 9' 'Y'Y' '9' 



,0), 



$^;C,/Z3 ■■ ^2, - • • , Sn) = (12; 0, 2, 4, 8, 14, 20, 28, 36), $3 .^3, and $1/2. 

m 12 4 4 4 4 4 4 ^s..Za '^^ old Zg parafermion state. $9.^(2) and 

( ^ 1 /^i > ' ^ ~ ^ 9 ' 3' 3' ^' 3' 3' 3' 3'^^' $_9_.^(4) are new Zg parafermion states, which can be 
(no, - - - , n„_i) = (2, 0, 2, 0, 1, 0, 2, 0, 2, 0, 0,0). 
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expressed as a correlation of simple-current operators ip2 
and ^4 in the Zg parafermion CFT, respectively. We 
also note that the /i-vector for the ^^-Cg/Za state is the 
difference of the h-vectors of the $ 9 state and the 
$3.7 state: 

4 4 44 44 
^3'3' '3'3' '3'3' ^ 

= (2, 2, 2, 2, 2, 2, 2, 2, 0) - (|| 0, II 0, II 0). 

However, we do not know if a symmetric polynomial de- 
scribed by the h vector (-^; |, |, 2, |, |, 2, |, |, 0) really 
exists or not. 



provide a fairly complete classification of non-chiral topo- 
logical orders in 2 spatial dimensions. We hope the 
framework introduced in this paper be a step towards 
a classification of chiral topological orders in 2 spatial 
dimensions. 
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IX. DISCUSSION 

In this paper, we use a local condition - the pattern of 
zeros to classify symmetric polynomials of infinity vari- 
ables. We find that symmetric polynomials of n-cluster 
form (see (23), (24)) can be labeled by a set in integers 
(m; ^2, • • • ,Sn). Those integers must satisfy the condi- 
tions (43), (45), and (47). 

Using the symmetric polynomials labelled by n and 
(m; • • • ,Sn),we have constructed a large class of sim- 
ple FQH states. The constructed FQH states contain 
both the Laughlin states and non-Abelian states, such 
as the Read-Rezayi parafermion states, the new gen- 
eralized parafermion states, and some other new non- 
Abelian states. Although, the constructed FQH states 
are for bosonic electrons, the bosonic FQH states and 
the fermionic FQH states have a simple one-to-one cor- 
respondence: 

^fermion — ^boson J_ J_ (^i ^i)* 

i<j 

We can easily obtain fermionic FQH states from the cor- 
responding bosonic ones. 

We have seen that the ground state wave functions of 
different Abelian and non-Abelian fraction quantum Hall 
states can be characterized by patterns of zeros {Sa}- 
One may wonder can we use the data {Sa} to calcu- 
late various topological properties of the corresponding 
fraction quantum Hall state? In Ref. 49, we will show 
that many topological properties can indeed be calcu- 
lated from {Sa}, such as the number of possible quasi- 
particle types and their quantum numbers. 

However, {Sa} cannot describe all FQH states. More 
complicated "multi-component" FQH states, such as = 
2/5 Abelian FQH state, are not included in our con- 
struction. This suggests that certain non-Abelian state, 
such as the parafermion states belong to the same class 
as the simple one-component Laughlin states. Thus, 
our result can be viewed as a classification of "one- 
component" FQH states, although the precise meaning 
of "one-component" remains to be clarified. String-net 
condensation and the associated tensor category theory^^ 



APPENDIX A: ELECTRONS ON A SPHERE 

To test if the FQH Hamiltonian has an energy gap or 
not, we need to put the FQH state on a compact space 
to avoid the gapless edge excitations which are always 
present. In this section, we will discuss how to put a 
FQH state on a sphere. We assume that there is a uni- 
form magnetic field on a sphere with a total iV^ fiux 
quanta. The wave function of one electron in the first 
Landau level has a form^^ 

^(61, ip)=Y, Cmv^i'-"'vr (Al) 

m=0 

where u = cos(6l/2)e'^/2 ^nd v = sin(6'/2)e-*'^/2.5i 

We can also use a complex number z = u/v to param- 
eterizes the points on the sphere. In term of z, the wave 
function becomes 

^{6, if) = e-'^*^/2sin^*(6'/2) ^ c„^™ 

m=0 

= (l + |^|2)^./2 l^ ^rnZ"" (A2) 
^ 11/ m=0 

We see that we can use a polynomial $(z) to describe the 
wave function of one electron in the first Landau level: 

^{z) = g c„^™ (A3) 

m=0 

Here the power of z is equal or less then N^. (Al) and 
(A2) allow us to go back and forth between the spinor 
representation 'if{u,v) and the polynomial representation 
$(^) of the states on sphere. 

Since the polynomials (A3) represent wave functions 
on a sphere, hence they form a representation of SU (2) 
(or 0(3)) rotation of the sphere. The dimension of the 
representation is Nff, + 1. Such a representation is said to 
carry an angular momentum 
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The SU{2) Lie algebra is generated by 

L'' = zd^-J, L-=d^, L+ = -z'^d^ + 2Jz, 
which satisfy 

[L^L=^] = L±, [L+,L-]=2L\ 

Those operators act within the space formed by the poly- 
nomials of the form (A3). The inner product in the space 
of the polynomials is defined through the inner product 

of the wave functions ^i{6,^p) 



77-72 $i(^) and 



($2|$i) = j sm{0)dedip %ie,ip)-9i{e,ip) 

-I 
-I 



4cos(^/2) sm^{e/2)d^^^^j^d<fi 



M^z 



sm{e/2) 



(1 + |2|2)2 (1 + |^|2)JV^ 



Now let us consider a polynomial of two variables 

^{^1,22) where the highest power for zi is 2Ji and the 
highest power for Z2 is 2J2 (here 2Ji and 2J2 arc in- 
tegers). $(zi,Z2) can be viewed as a representation of 
SU{2) too, where the generators of the SU{2) Lie alge- 
bra are given by 





= LI + L2, 


i± = 


-Lt- 


^Lt, 




LI 


= zidzi - Ji, 


Lr = 




L+ = -zfdz, - 


F 2JiZi 


LI 


= Z2dz2 - 






Lt = —Z29z2 - 


h 2J2Z2 



^(-21) -^2) is not an irreducible representation of SU{2). 
It can be decomposed as ©jt^jj^.j^i^j where is an 
angular-momentum- J representation of SU{2). We may 
say that z\ has angular momentum Ji and Z2 has angular 
momentum J2. So the angular momenta of $(zi, Z2) are 
those obtained by combining the angular momentum J\ 
and the angular momentum J2. 

What are the states in the space Hjl Let $j,m, fn = 
— J, — J + 1, • • • , J, be the polynomials in the Hj space 
such that is the eigenstate of with eigenvalue m. 
Let us also introduce z± = z\±Z2- We see that 



1 



2d 



L+ = -t{zl + zl)d,^ 



- z+z^dz_ + 2JiZi + 2J2Z2, 
Ji - J2. 



From L ^j-j = and L^^j^^j = -J^j^^j, we find 
that 



DC Z 



J1+J2-J 



{zi - Z2) 



J1+J2-J 



(A4) 



^j,m's are generated from $,7.-,/ by applying L+'s. Since 
never reduce the power of Z-, thus has a form 



zJi+J^-J f(zi,Z2) = (zi - Z2Y'+'^^--^ f{zi,Z2). This re- 
veals a close relation between the order of zeros and the 
angular momentum: the polynomial $(zi, Z2) with angu- 
lar momentum J must have an order J\ + J2 — J zero as 

Zi Z2. 

On the other hand, let ^d{z\, Z2) bo a polynomial that 
has a D^^ order zero as — > ^;2, «e has a form {z\ — 
Z2)^ f {zi , Z2) ■ We note that the actions of and 
do not decrease the power of zero as zi Z2. So the 
action of and can never change to ^j^-j if 
J > Ji + J2 — D. Thus ^d{zi,Z2) can only contain 
angular momenta J < J\ + J2 — D which lead to zeros of 
order D or more. 

We can let z\ Z2 in (-21,-22) and obtain i^(-22) as 



^d{zi,Z2) = {Zi - Z2)^-F(Z2) + 0[{ZI - Z2) 



The minimum eigenvalue for {z\ — Z2)^F{z2) is 
—J\ ~ J2 + D which corresponds to -F(.22) = 1- After the 
SU{2) rotations, we generate other polynomials F{z2) 
from F{z2) = 1- Those i^(.22) polynomials form an an- 
gular momentum Ji + J2 — D representation of SU{2). 
Thus fusing z\ and Z2 together produces a variable with 
an angular momentum Ji + J2 — D, where J\ is the an- 
gular momentum of zi , J2 is the angular momentum of 
Z2, and D is the power of the zeros as Zi Z2- 



APPENDIX B: THE RELATION TO 
CONFORM AL FIELD THEORY 

It was pointed out that the symmetric polynomial 
^{z\, - ■ ■ ,zn) that describes a FQH state can be written 
as an A^^-point correlation function of a certain operator 
Ve in a conformal field theory (CFT)!^'^!'^^ 

mzi})= hm z'^-{V{z^)l[V,{zi)) (Bl) 



To find the corresponding CFT of a symmetric polyno- 
mial <I> described by a pattern of zeros {Sa}, we will first 
calculate the scaling dimension of the operator Ve from 

Sa. 

1. The spin of the type-a particles 

Due to the relation between the scaling dimension and 
the intrinsic spin, we will first calculate the intrinsic spin 
of the electron and the a-electron clusters. From (8), we 
see that the angular momentum Ja of the type-a parti- 
cle has a form Ja = aJ — Sa, where Sa is an intrinsic 
property of the type-a particle and is independent of the 
magnetic flux A^^ = 2J through the sphere (sec (33)). Sa 
only depend on the pattern of the zeros and is called the 
orbital spin of the type-a particle. ^^■'^^'^^ 

As discussed in Ref. 47, the orbital spin contains two 
contributions Sa = S^ + ha- S^ comes from the spin 
vector and ha is the intrinsic spin. The intrinsic spin is 
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related to the statistics of the particle through the spin- 
statistics theorem. The statistics in turn is related to 
the scaling dimension (for example, bosons always have 
integral scaling dimensions). 

To separate the two contributions, we need to iden- 
tify the contribution from the spin vector. This can be 
achieved by noting that the spin vector contribution is 
proportional to o: = Ca. The key is to find the 
proportional coefficient C. 

For this purpose, let us consider SnN^ in (38). SnN^ is 
the orbital spin for the bound state of Nc type-n parti- 
cles. We know that the typc-r?, particles form a Laughlin 
state (20). For a bound state of type-n particles, its 
orbital spin SnN^ contains a term linear in A^c which is the 
contribution from the spin vector and a term quadratic 
in Nc which is the intrinsic spin. From (38), we see that 
the contribution from the spin vector to SnNc is 

TTiTi 

After replacing Nc by a/n, we identify the contribution 
from the spin vector to Sa- 



CSV _ „(^n '^x 



(B2) 



Thus the intrinsic spin is 



ha = S.-a{^-'^). (B3) 
ha is also the scaling dimension of the operator (K)"- 



The intrinsic spin ha is actually the scaling dimension 
of the a*'' power of the electron operator, Va = (14)"- 
The scaling dimension of the Gaussian part e""^^^^/^ is 



^-'t = ^and 



ha 



2n 



(B4) 



is the scaling dimension of the simple-current operator 

V'a- 

We can obtain the scaling dimension ha of the oper- 
ator Va more directly without using the concept of spin 
vector and orbital spin. First we note that the derived 
polynomial -P({^|"^}) can be expressed as a correlation 
ofK(4"^)'s 

P{{zt^})= hm z^^^{V{z^)\{Va{zt^)). 

The operator product expansion of VaS determines the 
powers of zeros in the correlation function P{{z["'^}). The 
{DabY^ order zero as 2;|"^ — > zf^ implies that 

Va{zi)Vb{z2) - (Zi - Z2)''-'Va+b + 0{{Zi - ^2)^- + ^) 

Thus the scaling dimension of Va satisfies 

ha+b = ha + hb + Dab (B5) 

From (10), we sec that ha and Sa satisfy the same equa- 
tion (except the = condition). Thus 



2. Symmetric polynomial as a correlation in a 
conformal field theory 

The electron operator Ve{z) in the GET expression of 
$ (Bl) has a form 



Ve{z) = Mz)e 



where e"^/^ is the vertex operator in a Gaussian model. 
The vertex operator has a scaling dimension i^. tjji is a 
simple-current operator, ie ipi satisfy the following 
fusion relation 

tpatpb = tpa+b, 'fpa = (V'l)°- 

Such an Abelian fusion rule is closely related to the 
unique- fusion condition discussed in section III A. If 
$(^1, ■ • • ,zn) has an n-cluster form, Vi satisfies 

Vn = (V'l)" ~ 1- 

<l?(2:i,-- - ,zn) can be decomposed according to (21). 
The correlation of the Gausian part e^"^^^^!^ produces 
the part of $ and the correlation of the simple-current 
part produces the G part of $. 



ha = Sa- Ca 

for a certain constant C. Since tpnk = (V'l)"'^ = 1 
and Vkn = e'"'^*/^, hnk contains only a term that is 
quadratic in k. From (38), we see that 



Snk — kSn + 



mnk{k — 1) 



(B6) 



Thus we must choose C to cancel the linear k term in 



(B6) to obtain hnk- We find that C = 



_ aSn am 
n 2 



and 



(B7) 



This allows us to obtain the scaling dimension of the 
simple current ijja which is given by h'^a in (34). We can 
also express Sa in terms of h^a (see (35)). 

We note that CFT requires that ha = mod 1, since 
Va are bosonic operators. This requires 



n 2 



(B8) 



Such a condition is satisfied if we impose the condition 
(51). 
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Let us introduce V-a = where we have assumed 

ipl = i^n-a- Consider 

G{zu...,Zk) = {Va,{^l)---Va,{zk)) 

i<j 

where Da.fc = ha+b — ha — hi,, P is a polynomial of Zi's, 
and X^jfl, = 0. The part ni<j('2^j ~ ^3)^°"'' reproduces 
the poles/zeros of the correlation function as Zi — Zj 0. 

As zi 00, Y[i<j{zi — Zj)^"''' behaves as z^ '"^ . As 
a CFT correlation function G{zi, . . . ,Zk) should behave 
as l/zi^"^ in 2:1 —> 00 limit. Thus, the maximal power of 
zi in P{zi, ■ ■ ■ ,Zk) must be 



7l = -^Da^^ai - '^hai, 
i=2 

k k 
= {k- i)hai + ^hai - ^ hai+ai, 
i=2 i=2 

where ai = — X]i=2 '^i- When fc = 4, we have 

Tl ~ ha2+a3+a4~^ha2~^ha3~\~ha4 ^a2+a3 ^a3+a4 ^a3+a2) 

where we have used h^a = ha- The requirement that 
7i > is the third condition in (45). 

3. Generalized vertex algebra 

To understand the CFT representation of the symmet- 
ric polynomial more deeply, let us consider generalized 
vertex algebra.^" The CFT formed by the simple-current 
operators ipa is a special case of a generalized vertex al- 
gebra. 

Consider operators A(z), B{w), etc which form an 
operator-product-expansion algebra 



A{z)B{w) = {\AB]^M 



(B9) 



+ {z- w)[AB]„^^_^{w) + {z- wf[AB]„^^_2{w) + • 
and 

{z - wy''A{z)B{w) = hab{w - z)"'"'B{w)A{z), 

(BIO) 

where /kab is a phase factor. Here {z — w)"-*^ = \z — 
^^aAB (.aABiO^^ where (z — w) = \z ~ wje'^^™ and — tt < 
0ZW < TT. The operator product in (B9) is assumed to be 
radial ordered: A{z)B{w) — > R[A{z)B{w)] where 

iz-w)"^''R[A{z)B{w)] 

_ j{z- w)''^''A{z)B{w), \z\ > \w\, 

~ \ fJ-Asiw - z)"'"'B{w)A{z), \w\ > \z\. 



We see that the commutation relation (BIO) ensures that 
the correlation functions of A{z) and B{'w) are smooth 
functions. Let /ia, hs and h[AB]c^g be the scaling di- 
mensions of A, B, and [ABj^j^^, then 

CXaB = hA + hs — h\^AB]aj^jj ■ 

The self consistancy of the vertex algebra requires 
a^s's to statisfy^^ 

ttAB + CKyic - "AD = mod 1, (Bll) 

and fiAB^s to statisfy 

t-CABt-CAC = /XAc(-)"^"+"^^-"^° , (B12) 

where D = [BC]aBc- 

The fusion rule of the simple-current operartors 
V'aV'6 = V'a-i-6 requires that those operators form the fol- 
lowing vertex operator algebra: 



%pa{z)lp-a{w) 



Cab{^a+b{w) + 0{Z - wfj 

f \/iSC_l_/,SC_l,SC 

[z - wy" 



1 + ^{z - wfT{w) + 0{{z - wf) 



{z — w) 



2h%' 



where tp-a = tpn-a and tpa+n = tpa- 

We see that the above algebra of simple currents tpa is 
a special case of generalized vertex algebra. We have 

aab = h^a + hT - h^a+b- 
The condition (Bll) becomes 

A3(o,6,c) =0 mod 1, (B13) 

where we have used (49) and (46). From (B12), we see 
that if A3 (a, b, c) = odd for certain choices of a, b, c, then 
fj,ab cannot be trivial {ie jiab = !)• When A3(a, 6, c) = 
even for all a, b, c, then ^ab can be a trivial solution fXab = 
1. Thus the condition (51) has a special meaning in CFT. 



4. The conditions on the /i-vectors 

Due to the one-to-one conrespondcnce between the S- 
vectors and the vectors [see (34) and (35)], we can 
translate the conditions (43) and (45) on the S-vectors 
to some conditions on the /i-vectors. 

Note that a /i-vector is specified by n, m, and 
/if , • • ■ , h^-i- We extend /if , • • • , /i^^_i to h^a for any 
integer a by requiring 



= 0, 



"a "'a+n- 



The conditions (43) become 

ba = — ahi -\ = non-negative mteger. 



2n 



m > 0, mn = even, 
2n/if -I- ro = mod n. 



(B14) 
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From 2nh\^ + m = mod n, wc see that 2nh^^ is an 
integer. From 2nh^^ —a{2nKf )+a{a—l)m = even integer, 
we see that 2nh^ are always integers. Also 2n/i2o 
always even integers and 2n/i|^_,_j are either all even or 
all odd. Since = 0, thus when n = odd, 2nh^ are 
all even. Only when n = even, can 2n/i2o+i either be all 
even or all odd. 

To generate sets of that satisfy the above condi- 
tions, we will use (34). Setting a = 1 in (34), wc get 
2nhf' = m{n — 1) — 2S'„. We see that when n = odd, 
2nh^^ = even. When n = even, 2n/if'^ = even when m = 
even and 2nh^^ = odd when m = odd. We also see that 
Sn < m{n — l)/2 which implies that. Sa < m{n — l)/2 
for a = 2, 3, • • • ,n. 

The conditions (45) become 



2a 



LSC _ tsc _ ti 
"■a+5 " 



2anhf + ma{2a - 1) 

^ ahm 
b ^ 



mod 2n, 



> -- 



LSC _ LSC _ LSC _ tSC i ts< 

"'a+b+c "'a+b "'b+c "'a+c ' "'a 



+ hi'' + K 



(B15) 
> 0. 



The condition (51) becomes 



K+b+c - K+b - K+c - K+c + K + K + K = even. 

(B16) 



The first condition in (B14) and the second condition 
in (B15) imply that 



^SC _ LSC _ 

"a+6 "a 



n 
ahm 



K+b - K -hf + — = Q mod 1, 



(B17) 



which are part of defining conditions of parafermion CFT 
if m = 2. So the pattern of zeros of symmetric polynomial 
may have a natural relation to parafermion CFT. 
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